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ABSTRACT 


This  report  describes  a  modified  ray-tracing  technique 
used  in  the  synthesis  of  oblique -incidence ,  step-frequency 
ionograms.  Ionograms  of  this  type  are  obtained  experi¬ 
mentally  to  aid  in  the  real-time  selection  of  frequencies 
for  point-to-point  communications  and  propagation  studies. 
When  it  is  desirable  to  identify  the  modes  of  propagation, 
computer-calculated  ray  tracings  have  proved  quite  valuable. 

The  Kift-Fooks  ray-tracing  technique  was  chosen  because 
it  is  a  rapid  program  capable  of  tracing  rays  when  only  a 
minimum  of  ionospheric  data  is  available.  One  could  utilize 
this  technique  in  the  analysis  of  propagation  data  either 
by  synthesizing  an  oblique-incidence  ionogram  for  direct 
comparison  with  experimentally  observed  results  or  by  com¬ 
paring  plots  of  maximum  usable  frequency  (predicted)  with 
receiving-station  log  sheets.  The  details  of  the  computer 
program  are  included  with  instructions  that  may  be  used  as 
a  guide  by  anyone  familiar  with  computers  and  programming 
operations  to  perform  his  own  calculations . 
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t,  INTRODUCTION 


In  an  effort  to  understand  better  the  propagation  char¬ 
acteristics  associated  with  fixed-frequency  transmissions 
over  a  long  (8000-km),  cast-west  path  (l.e.,  Hawaii  to 
Massachusetts),  it  was  decided  to  instrument  this  path  with 
a  3 t.e p - f re que nc y  (4-64-Mc)  transmitter  and  a  synchronized 
receiver.  With  the  above  equipment  operating  on  a  round- 
the-clock  basis,  it  was  hoped  that  records  could  be  obtained 
that  would  permit  deduction  of  the  mode  structure  and  appar¬ 
ent  ray  path  of  the  propagating  signals. 

Examination  of  'he  records  taken  on  this  path,  soon 
indicates  that  the  usual  simplifying  assumptions  (such  as  a 
uniform  ionosphere  ov< r  the  entire  path)  are  often  not 
representative  of  what  is  happening.  The  path,  3000  km  long, 
is  Just  on  the  edge  of  the  normally  assumed  "allowable”  two- 
hop,  F2-layer  propagation.  Records  show  that  the  2F2  mode 
propagates  for  only  short  periods  around  noon  and  midnight, 
local  time,  at  the  midpoint  of  the  path.  At  other  times 
(especially  sunrise  and  sunset),  the  progressive  change 
across  the  path  from  a  daytime  ionosphere  (with  E,  Fl ,  and 
F?  layers)  to  a  nighttime  ionosphere  (F2  layer  only)  produces 
a  bewildering  variety  of  propagation  modes.  Analysis  soon 
becomes  fairly  complex.  To  assist  in  the  understanding  of 
the  mod**  structure,  it  was  felt  that  a  ray-tracing  program 
that  simulated  the  experimental  data  would  help. 

Familiarity  with  the  experimental  technique  and  record 
form  will  help  in  under standing  the  type  of  information  that 
would  be  desirable  from  Computed  ray  tracings.  The  experi¬ 
mental  data  were  obtained  in  the  following  way.  The  trans¬ 
mitter  and  receiver  include  electronically  tuned  and  syn¬ 
chronized  circuitry  that  ranges  in  frequency  from  4  to  64  Me? 
in  160  st^ps*4o  linearly  spaced  steps  per  octave  band. 

Pulses,  50*sec  in  deration,  are  transmitted  over  the  Hawaii - 

.  I  - 


Massachusetts  path,  and  the  received  pulses  (differentially 
delayed  In  time  as  a  result  of  the  different  modes  of 
propagation)  are  recorded  on  film  using  the  following  tech¬ 
nique,  An  oscilloscope  Is  Intensity  modulated  with  the 
detected  video  output  of  the  receiver.  As  the  transmitter 
and  receiver  step  in  frequency  over  the  operating  range,  the 
display  is  recorded  on  moving  film,  producing  a  record  show¬ 
ing  time  delay  as  a  function  of  frequency.  This  presentation 
is  referred  to  as  an  oblique -Incidence  lonograin.  An  artist's 
sketch  of  this  type  of  record  is  shown  in  Fig.  1. 

The  primary  characteristics  that  one  would  hope  to 
obtain  from  a  ray-tracing  analysis  for  comparison  with  the 
experimental  results  are  summarized  below: 

1.  The  Maximum  Observed  Frequency  (MOF)  and  the  Lowest 
Observed  Frequency  (LOF)  for  each  of  the  modes 
(e.g. ,  mode  1,2,3, . 

2.  The  differential  group  time  delay  separating  each 
of  the  modes  at  any  given  frequency  (e.g.,  f^). 

In  addition.  It  would  be  desirable  to  obtain  a  profile  viev: 
of  the  propagation  path  showing  the  rays,  their  ground- 
reflection  points  and  the  apparent  path  of  the  rays  through 
the  Ionosphere.  An  -  /ample  of  this  is  given  in  Fig.  2, 
showing  the  three  modes  of  the  ionogram  of  Fig.  1,  at  a 
fixed  frequency  f ^ . 
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II.  AVAILABLE  RAY -TRACING  TECHNIQUES 


The  problem  in  synthesizing  an  obi tque-incidence 
ionogram  by  a  ray-tracing  approach  i3  actually  twofold: 

1.  Can  the  mode  structure  be  duplicated  by  a  ray-tracing 
approach  if  sufficient  ionospheric  data  are  available? 

2.  In  the  absence  of  this  ionospheric  data,  could  the  CRPL 
ionospheric -propagation  predictions,  available  three  months 
in  advance,  be  used  in  conjunction  with  the  ray-tracing 
program  to  predict  the  mode  structure  likely  to  be  observed? 

With  this  problem  in  mind,  it  was  decided  first  to  find 
out  how  other  researchers  had  solved  this  or  similar  problems. 
Inquiry  into  the  available  ray -tracing  techniques  necessi¬ 
tated  visiting  various  establishments  to  find  out  the  latest 
information;  at  that  time,  much  of  it  was  as  yet  unpublished. 
However,  since  then,  a  meeting  has  been  held  in  Lindau, 
Germany,  to  discuss  oblique -incidence  soundings  and  iono¬ 
spheric  ray  tracing. 

Table  1*  13  a  summary  of  ray-tracing  techniques. 

An  alternate  possibility,  the  use  of  an  analog  computer 
to  solve  the  ray-tracing  equations,  ha3  boon  utilized  hy 
Wong  [Ref.  17].  The  difficulty  in  using  an  analog  computer 
la  that  the  output,  height  vs  range  (as  a  function  of  fre¬ 
quency),  gives  the  distribution  of  energy  along  the  great 
circle  but  doe3  not  "home -in"  on  the  receiver  (a  point  at  a 
fixed  range). 


This  information  is  based  on  material  that  appeared  in  the 
“Report  oT  the  Lindau  Meeting  on  Oblique  Sounding  of  the 
Ionosphere /  May  6-10,  1963.  Meeting  held  at;  Institut 
?iir  .lonsopharen-phyaik,  Max-Planck-Irstitut  Fur  Aeronomie  , 
Lindau  Nartheim,  Germany. 
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'the  !'home -iri“  capability  afforded  by  a  digital -computer 
program  Is  fen  attractive  feature.  Sorting  by  modes  and  range 
discrimination  greatly  simplifies  the  handling  of  the  enormous 
amounts  of  data  that  are  calculated  by  the  computer .  One  is 
thus  able  to  concentrate  on  the  path  In  question,  having  al¬ 
ready  sorted  out  the  rays  that  never  reach  the  receiver. 


fl<i.  \.  tlitORtTlcm  Of  4*  OBLIQUE'  INCIDENCE.  SW1UP 

fttfuttfcv  10NQGB4*. 
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III.  CHOICE  OF  THE  KIFT-FOOKS  TECHNIQUE 


To  synthesize?  an  oblique-incidence  i  one  gram  (Fig.  1). 
it  is  necessary  to  consider  only  thos-  rays  that  reach  the 
receiver.  Detailed  knowledge  of  the  ionosphere  is  not 
always  available  and,  wh-re  predictions  are  cons  rood,  a 
detail'd  ray-tracing  approach  is  not  justified,  la  fact, 
most  of  the  tltm  ,  only  a  car**  minimum  of  data  xisti  con¬ 
cerning  the  true  <  1<  ctron-u*  :.si ly  prof. 1.1 f  along  any  glV1  n 
path.  Ever,  with  el'ctron-d  nsit;,  distributions,  assump¬ 
tions  a.;  1  o  the-  structur  of  the  magnetic  fl*  id,  the  off- 
gr'-at-circl*  pr  files,  as  w  *11  as  a  choicf  of  a  magneto- 
ionic  th  ory ,  need  to  t  <  ra  h*  prior  to  th*1  use  of  a  'OR.pl-  t' 
three-dimensional  analysis  [R'f.  lr  ]  . 

With  the  a  limltat  ions,  .it  was  believed  that  a  program 
which  tak  s  Into  account  tic.  gross  changes  in  the  ionosphere 
along  a  path  at  sunrise  cu.J  sun::**t,  by  the  inclusion  of  the 
daytime  E  and  FI  layers  a..d  a  specularly  reflecting  sporadic 
E  layer,  viouj d  suffice. 

The  major  factors  gov-  ruing  the  choice  of  the  Kift-Fooks 
technique  were  probat ly  the  rapidity  with  which  the  program 
could  be  run  on  a  truly  high-speed  digital  computer  (oitn  r 
the  IBM7090  or  the  IRM  rOn4)  and  the  fact  that  predict  ion., 
could  be  made,  using  the  CRFL  ionospheric  propagation- 
predictions  in  their  pres-  ut  card  format  [R*  f,  lqJ,  on  a 
highly  automated  basis. 

Thus,  it  was  decided  to  use  the  ray-tracing  technique 
suggested  by  Kift  [Ref.  11]  and  programmed  for  use  on  the 
?egasus  computer  by  Fooks  [Ref.  12].  The  advantages  of  this 
program  are  that  It  assumes  a  set  of  parabolic  layers  for 
the  ionospheric  profile  and  then  calculates  the  ray  path  in 
(or  through)  a  parabolic  layet  by  the  Apple  to 
*  [Ref.  equations. 
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Some  of  the  inaccuracies  of  this  tachriiQue  ai'1  pointed 
out  by  Kift  at  the  end  of  !>is  article,  vith  reference  to  the 
work  of  Vickers  [Ref.  19J.  A  report  that  compares  the  Kift- 
Fooks  technique  with  a  more  accurate  technique,  developed  by 
Croft  [Ref.  20],  is  soon  to  be  published  as  another  r°port 

in  this  series. 


flG.  2,  CROSS  SPCflWN  Of  IONOSPHERIC  H^T  PATHS, 


*>  ^  * 
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IV.  HOW  TO  UT»IM2S  ¥HB  KIFT-FOORS  VECflNIQ.Ug 


Before  attempting  to  utilize  fcho  Xlfl-Foohs  technique 
in  the  analysis  of  point— Co— point  propagation  charaot^rlstlos, 
one  should  know  what  data  ar»  available  In  the  prlPfc-ott  of 
the  program  and  how  to  use  these  data. 

Table  2  lists  the  Ionospheric-profile  parameters  along 
the  great-circle  path  (from  the  transmittw  r  to  the  !*•  coiver) 
in  100 -km  intervals.  The  values  given  are  the  critical  fre¬ 
quencies  of  the  E,  PI,  and  F2  layers,  and  the  height  at  which 
th'  maximum  of  th*.  F2  layer  occurs  is  given  lor  tin  points 
mentioned  above.  Details  of  the  exact  computer  output  format 
(Tables  2,3,'*)  are  given  on  pp.  11,  12,  and  13. 

Each  ray -tracing  group  is  identified  by  the  transmitter 
latitude  and  longitude,  the  bearing  to  the  receiver.  Aid 
the  tin  ,  month,  and  year  for  which  the  ionosphere  was 
complied . 

An  exanpl  of  the  pr<  i<  nt  data  format  is  given  in  Tal  le 
3.  fiv-  description  of *'  ach  mod  includes:  names  ol  succes¬ 
sive  reflecting  layers,  frequency,  take-off  angle,  gro^p 
»  ir  -  rj.  lay,  and  aUwnuatio-  .  'hi  V  rminology  use- 1  in  this 
format  is  dlflVren,  fro;;,  that  r  comm-..iei  for  use  in  oell/U-e- 
incidenc"  woi'k  (Appendix  i  .  Howev  r,  since*  this  report  i.> 
intended  to  explain  th*  ray-tracing  program  in  its  present 
form,  inclusion  of  the  recommended  u.ot/  nclature  would  have 

necessitated  further  delays. 

The  modes  are*  listed  In  terms  of  ineroacln*  frequency 

and  take-off  angh  (for  any  one  frequency). 

In  addition,  an  option  available  to  the  program  pPir.ia 
out  the  ground  range  and  height  of  the  ray  for  point D 
entry  or  exit  of  a  layer  and  the  ground-reflection  points 
(Table  4).  Thus  a  ray  plot  similar  to  that  shown  In  Figure 
2  could  be  plotted  from  the  data  of  Tati' 

With  knowledge  of  the  output  format  in  hand,  on^  can 
now  proceed  U»  diceuawinn  of  how  data  ear. 

.  4?  .  SE&-6>!&3 


used  in  the  synthesis  of  a re  oblique -Incidence  ionogram 
(Fig.  l).  Re ff erring  to  Tabl1?  3  ah&  establishing  thSf  sairi- 
set  of  coordinates  as  that  achieved  experimentally,  one 
would  then  plot  and  join  together  points  having  the  same 
mode  description  (i.e.,  .FI  .5  ,E  .E  .E).  This  plot  eou_^ 
then  be  compared  directly  with  the  experimentally  achieved 
data.  Please  note,  however,  that  th-  re  Will  he  a..  *»•  -- -- 
of  the  high -angle  rays  because  of  the  method  use  d  pro- 

grammir.g  the  computer  for  mode  calculation  and  retention. 

When  the  take-off  angle  and  the  attenuation  associated 
v.’ith  a  given  mode  are  taken  into  account,  a  i^rs^-order 
approximation  can  be  made  to  eliminate  many  of  the  predicted 
modes  that  experience  tells  us  just  wouldn't  get  tnrough. 

Using  vertical-incidence  soundings  made  along  or  near 
the  great  circle,  as  a  first-order  correction  to  the  CRPL 
prediction,  enables  greater  accuracy  to  be  achieved,  par¬ 
ticularly  if  patches  of  sporadic  E  are  present  which  were 
not  taken  into  ae  -o.nt  in  the  predictions.  A  subsequent 
report  will  be  issued  outlining  the  procedure  used  in  this 
case  (I.e.,  an  after-the-fact  analysis). 

However,  it  is  most  important  to  emphasize  once  again 
the  main  advantage-  „f  the  Klft-Fooks  tecnr;xOu‘.  a-  a  pr  - 
dictor  of  propagation  conditions.  Certainly,  when  detailed 
information  regariir.g  the  ionospheric  profile  is  available  , 
it  would  make  sense  to  utilize  one  of  the  more  detailed 
ray-tracing  programs  currently  available  [Refs.  ),l4,lo,*  . 

By  directly  converting  the  CRPL  ionospheric  propagation 
predictions  into  values  of  fQF2  and  *3003  (the  ratio  of  the 
3000  Km  MUF  to  the  fQF2)  and  subsequently  using  the  assump¬ 
tion  of  Kift  and  Fooks  [Refs.  11  and  12,  respectively},  tr 
computer  can  calculate  the  values  of  height  of  the*  maxifiiua 
of  the  F2  layer  and  trace  all  subsequent  rays  that  reach  he 

receiver1 . 

Thus  we  have  a  highly  automated  prediction  program,  the 

avails  of  wtii.cn  ate  ap'.UfUh  ...  tM  foliouir*  section. 
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TABLE  2.  IONOSPHERIC- PROFILE  PARAMETERS 


u 
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TABLK  3.  OBLIQUE* IONOGRAM  OUIPUT  DATA 


4  OCTOBER  1942  1TJ7.  3t  CHI  RRMOR/8EOFURO  R1TH 


PATH  1 

length 

80*5. 

35  AH 

IX 

IAT 

19.8 

"00fc 

•  f 

.9 

.E 

.F  2 

.E 

.f 

a 

a 

a 

,C 

.E 

.E 

.[ 

•  E 

•  E 

.€ 

a 

.9  2 

.E 

a 

a 

a 

a 

.f 

a 

•f 

a 

.E 

•  E 

a 

.1 

•  t 

•  E 

a 

.f  2 

.92 

.92 

.92 

.92 

.92 

.92 

•  1 

•E 

•  1 

•  t 

.9* 

•»l 

•  *X 

•M 

.1 1 

•92 

•  92 

•92 

.9  2 

.92 

.12 

.12 

•92 

•92 

•  12 

•  92 

•  *2 

.9  2 

. 9 ? 

•  9 1 

•( 

•  1 

.9 

•  9 

.9  2 

•92 

•  92 

.92 

•  92 

.12 

•92 

•92 

•  92 

•*2 

•  H 

.12 

.12 

•92 

•92 

•  92 

•92 

•  12 

.12 

.12 

•92 

•  *2 

.91 

•  92 

•12 

•  >2 

•  92 

•9  2 

.92 

•U 

.12 

•  92 

•92 

•  92 

.92 

•  12 

.12 

.12 

•9  2 

•92 

•92 

•*2 

.*  2 

.12 

.12 

•  ' 2 

•*2 

.1  2 

•  n 

•  12 

•92 

•92 

•92 

.92 

•  *i 

•12 

•92 

• 92 

•  92 

•  92 

•  U 

.12 

.12 

•  92 

.92 

•92 

.92 

•92 

.92 

•92 

•  92 

.•2 

•92 

•92 

•92 

•  92 

•  12 

.12 

•  92 

•92 

•  92 

•  *2 

•*2 

.12 

.12 

•  92 

•92 

•  92 

•  *2 

.92 

•92 

.91 

•  92 

•  12 

•  92 

•  92 

•  92 

.9  2 

•*: 

a: 

•92 

.92 

•  12 

•  *  2 

.92 

.9  2 

.9  2 

.92 

.92 

•12 

.92 

.92 

.9  2 

.92 

.92 

.9  2 

.9  2 

.9  2 

.92 

.92 

.9  2 

.9  2 

.92 

.9  2 

.92 

.9  2 

.92 

.92 

.9  2 

.92 

.92 

.92 

.92 

.9  2 

.92 

.92 

.92 

.92 

.92 

.92 

.9  2 

.92 

.92 

.92 

.92 

.92 

.92 

.92 

.92 

.92 

.92 

.92 

.92 

.92 

.9  2 

.92 

.92 

.92 

.92 

.92 

T*  LONG  -1.4.9B  Ot.  M  BERR1B*  %'  .ft  If* 

IREO  BC«A  ulRI  t(“l  It*  » 


«. 30 

1.18 

8039.09  ft. 10 

•t m2%  649. 8? 

4.00 

1  3.f  9 

8041.41 

27.96 

•V>4  #9%.  86 

8.00 

1.22 

8037.90 

27.09 

•7.49 

r  .1.2  a 

6.X 

f  •  26 

8*136.66 

27.09 

•8.69 

hM* 

7.00 

1.35 

0038.48 

27.09 

»9.af  2H.a» 

7.00 

12.80 

8026.50 

2  7*84 

*18.89  2?8.9| 

0.00 

1.44 

8034. 8| 

2i.09 

•10.78 

163.94 

9. CO 

1.54 

8036.3* 

27.13 

•11.01 

129.5* 

10.00 

1.68 

6038.4) 

2  7.12 

•  9.92 

r*.n 

10.00 

2  3.19 

801  5.7* 

30.17 

-29.59 

*6.37 

li.oo 

1.88 

00)9.26 

27.16 

-6.C9 

00.  *6 

il.OO 

13.20 

7989.22 

26.10 

-58.13 

t'.J* 

11.00 

21.7C 

6024.16 

29.07 

-21.17 

•1.77 

11. GC 

24.17 

K41.8* 

30.51 

-1.51 

IM1 

12.00 

7.01 

8043.01 

27.81 

-2.12 

71.66 

12. Of 

18.43 

0021.69 

29.25 

-23.6% 
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V. 


© 

THE  KIFT - FQOKS  RAY  ^RhClhG_tROGBAK 


The  ionospheric  ray-tracing  program  describe d  h<  re 
essentially  the  same  as  that  described  by  G.  F.  Fooks  in  his 
r(  port  [Ref.  12].  The  same  equations  are  used  and  the*  name 

basic  procedure  is  followed:  however,  certain  oodlflcati  ns 
and  additions  to  the  program  have  been  made  to  allow  t.  e 
calculations  of  reflection  heights  from  the  Ionospheric 
lasers  and  to  allow  the  calculation  of  ar.  approximate  val 
for  x*ay  attenuation  due  to  D -layer  absorption  along  the  paU  . 


A.  PHYSICAL  ASSUMPTION.! 

The  program  uses  a  curv  J -earth,  curv  i-ionosph~re 
geometry,  and  the  ionosphere  is  assumed  to  consist  of  a 
number  of  curved  layers,  each  with  a  paratolic  electron- 
density  distribution.  The  ionospheric  layers  considered  are 
H  F  ana  F„.  A  sporadic  E  layer  (E=)  may  also  be  Ir.c  u  e 
lr.  the  calculations:  however,  when  It  is,  it  Is  treated  not 
as  a  oarabolic  layer,  but  rather  as  a  thin,  specular  ,  re. 
ting  ‘sheet.  The  earth's  magnetic  field  and  layer  tilts  are 

ignored • 

■  Figure  3  111  .strafs  the  geometry  of  the  Ionospheric 


layer  structure. 


f 

FIG.  3.  IONOSPHERIC- LAYER 
SIRUCi'L'RE  iPARAROl.TCio 
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R.  GENERATION  OF  THE  IONOSPHERE 


Figure  A  illustrates  the  geometric  parameters  foi  an 
arbitrary  parabolic  layer. 


110.  1.  OKOMHIUC  I’AltAMKIKItS 
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V  =  critical  t'r  queney  For*  th  layi'  in  iivga cycles 
(vertical,  incidence) 


SON  =  cun spot  number 
^  *  ffolar  angle. 
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For  the  F  layer: 


f  FI  =  1.4(f  E) 
o  o 

h  F  =  210.0  km) 
m  1 

y  F.,  =  60.0  km 

1  J 


Por  the  F^  layer,  values  of  fQF2  and  are  supplied  to 

the  program  either  as  predicted  values  or  observed  values 
at  arbitrary  points  along  the  path,  and  the  program  con- 
tructs  a  parabolic  F_  layer  under  the  assumption: 


y  F  =  0.4  h  F  , 
2  o  2 


(3) 


where  ym  =  t>m  -  h0 

Values  of  f  Eg,  if  they  are  different  from  zero,  are 
supplied  to  the°program  in  terms  of  their  position  on  the 
path.  The  height  of  the  Es  layer  is  assumed  constant  at 

100.0  km. 

An  equation  for  cos  x  using  the  path  geometry  is  pre¬ 
sented  In  Appendix  B. 

In  Appendix  C  a  method  is  given  for  obtaining  values  of 
h  F2  using  predicted  values  of  fQF2  and  F2  4000  MUF.  These 

are  the  two  parameters  obtained  from  the  CRFL  ionospheric 
predictions . 
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c.  EQUATIONS  FOK  RAY-PATH  CALCULATIONS 


Below  the  ionosphere  and  between  ionospheric  layers  the 
ray  is  assumed  La  travel  In  a  straight#  line. 

The  Ss  layer  either*  specularly  reflects  the  ray  or  allows 
it*  to  pass  dndevlated.  For  the  parabolic  layers  the  follow¬ 
ing  equations  apply: 

2f  <* 

AP'  =  Y~  'vm  *  ar6ta.nh  (jr  cos  i. )  (A) 

°  o 

if  the  ray  is  reflected  by  the  layer,  and 


2f  f 

aF’  =  r“yrn  ar6c°fh  (4  COG  i) 


o 


(5) 


if  the  lay  r  transmits  the  ray,  but  eaus'-o  bending,  where : 


AP 1  —  virtual  path  in  the  layer 
f  =  wave  frequency 

=  laYer  critical  frequency 


line  to  tb; 
tin  vertical 


■  i  t-he  ray,  extrapolated  along  a  straigh 
lev-.l  of  maximum  el.  ctron  density,  and 
at  that  leve  l  (as  illustrated  in  Fig. 


). 


riu,  5.  OiSl.lQIin- I  NCI  l)KN(  !•  HAY -PATH  GKOMKTRY. 
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For  transmission  through  a  lays?  ox*  *ox  £*efl5ict*ioft  frOT' 
the  bottom  of  a  layer*. 


AD 


■  R-TTT  slrl  1  AF‘* 

m 


If  the  ray  is  reflected  from  the  top  of  a  layer # 


AD  =  R  j  sin  1  A?*, 

m 


16) 


(7) 


where  AD  is  the  range  along  the  path  covered  while  the  ray 
is  in  the  layer,  and  R  Is  the  earth's  radius. 

In  the  course  of  the  ray  tracing,  as  the  ray  enters  a 
layer,  there  are  three  possible  consequences : 

1.  The  ray  is  reflected  from  the  layer. 

2.  The  ray  is  transmitted  through  the  layer  and  deviated. 

3.  The  ray  is  transmitted  through  the  layer  ur.de viated 


(straight-line  transmission). 


Let 


Then  if 

K  <  1 

K  =  1  F 1  =  00  > 

1  <  K  <  2 
K  1  2 


K  =  (f/fQ)  .  cos  i. 


the  ray  is  reflected 
the  next  ray  is  taken 
the  ray  Is  transmitted  and  i-  viated 
uh*.  ray  Is  brans:  it  ted  ar.J  undeviat-  d j 
The  equations  use  l  for  utri-  v latr  1  transmission  through 
a  layer,  between  layers,  and  free,  the  ground  the  uottor 
of  the  ionosphere  ar  : 


>  (3) 


Sih  1, 


( r  +  )  *-!.»♦  i- 


(9) 


(R  +  h,) 


AP 


’R  +  hr)  sin  ( i^  -  i?) 


X1 


AD  = 


_  R(in  -  io) 


(10) 


ClD 


for  straight-line  transmission  between  two  points  at  heights 

h  and  hQ  fts« opiated  vertical  anglsa  i-,*  a'«d.  i 

12  **» 


2* 
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During  the  course  of'  the  ray-tracing  procedure,  two 


layers  may  haopen  to  overlap  (most  likely  the  FI  and  F2 
layers).  When  this  occurs,  the  ray  is  extrapolated  back 
along  a  straight-line  path,  tangential  to  its  direction 


when  it  emerges  from  the  first  layer, 


to  its  point  of  entry 


to  the  second  layer,  (Fig.  6). 


LAYER  2 


LAYER  I 


HO.  <>.  0\ Kill  - AIM’ I N(i» l  AY Kll  INIOCKWHK 
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VI. 


STANFORD  VERSION  OF  KIFT-FOOKS  RAY -TRACING  PRO GRAF 


A.  BASIC  COMPUTATIONAL  PROCEDURE 

Figura  7  Is  a  logical  flow  diagram  of  the  computational 
procedure;  it  is  not  intended  as  a  detailed  flow  chart  oi 
the  program,  but  merely  as  a  gross  logical  description  of 

the  computational  process. 

The  input  data  to  the  program  are  the  path  length 
between  the  receiver  and  the  transmitter,  the  coordinates  of 
the  transmitter,  the  true  bearing  of  the  receiver  from  the 
transmitter,  sunspot  number,  sun  declination,  apparent  solar 
time  at  Greenwich;  F. -layer  data  lr.  the  form  of  either  fQF£ 
and  h  F  or  f  F,  and  ?,,  WOO  KUF ;  E,  data,  if  any,  plus  a 

range  of  frequencies  ar.d  a  range  of  take-off  angles  to  be 
Investigated  Tor  the  giver,  ionosphere,  and  a  set  of  fre¬ 
quencies  for  which  ray -reflection-height  information  is 

desired.  F-  and  data  are  described  lr:  terms  of  their 
d  ° 

range  along  the  path  from;  the  transmitt'  r. 

Once  the  data  for  the  path  have  been  read  by  the  program, 
a  table  of  ionospheric  uata  is  produced  lor  ust-  tl  th> 
program.  Equations  (l)  and  (2)  ax-  evaluated  at  If  -km 
intervals  along  the  path.  A  second- degre.  polynomial  is^ 
fitted  to  successive  triplets  of  iQF2  arid  h^F^  data  points 

and  these  polynomials  ar  evaluated  at  100-km  interv^.s 
along  the  path.  Tables  of  Eg,  if  required,  are  complied  at 
10-km  intervals  within  each  E^  patch  considered.  When, 
during  ray  tracing,  values  within  the  ionospheric  tables 
are  required  between  the  calculated  100-km  (10-km  for  Es) 
points,  linear  interpolation  is  used. 

After  the  generation  of  the  ionospheric  tables,  the  values 
of  the  critical  frequency,  fQ  and  height  vs  range  for  each 
layer  are  printed  ©ut  LTable  2)  and  the  actual  I,ay-tr*acin6 
process  begins# 

-  20  - 
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FIG.  7  B  fl.OW  lit  AG  RAM  OF  RAY» TRACE  PROGRAM. 
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For  each  frequency  of  the  specified  frequency  range, 
rays  are  traced  from  one  transmitter  for  all  take-off  angles 
of  the  specified  take-off  angle  range.  A  ray  is  traced  until 
it  falls  within  ±1000  km  of  the  receiver,  each  reflection 
from  a  layer  being  recorded,  in  coded  form,  in  a  "reflection 
index".  A  maximum  of  ten  reflections  per  ray  is  allowed. 

When  a  ray  falls  beyond  the  receiver  range  +1000  km,  tracing 
of  that  ray  is  terminated,  the  reflection  index  and  accumu¬ 
lated  values  of  F '  and  D  are  restored  to  zero  and  the  next 
ray  of  the  series  is  traced.  When  the  ray  falls  within 
±1000  km  of  the  receiver  range  its  reflection  index  is  stored 
in  a  table,  along  with  ?’  and  D  for  that  particular  ray,  and 
the  next  ray  of  the  series  is  traced. 

Once  all  rays  for  a  given  frequency  have  been  traced, 
the  table  of  reflection  indices  is  searched  for  rays  of  like 
modes,  a  linear  Interpolation  of  take-off  angle  with  respect 
to  the  actual  range  of  the  ray  and  the  range  of  the  receiver 
is  performed,  ar.d  a  new  ray  with  the  interpolated  take-off 
angle  is  traced.  At  this  time  the  height  of  the  ray  at  its 
reflection  points  is  calculated  in  addition  to  an  estimate 
of  D- layer  attenuation.  (These  equations  appear  in  Appen¬ 
dixes  D  ar.d  E,  respectively  .)  If  this  new  ray  does  not  fall 
within  ±100  km  of  the  r  ceiver,  it  is  ignored  and  the  next 
pair  of  like  rays  (If  any)  is  considered.  If  it  falls  with¬ 
in  ±100  km  of  the  receiver,  its  delay  time  is  calculated 


from 


Delay  Urn,  (...)  -  (12, 


and  the  final  results  are  printed.  The  next  pair  of  like 
rays  (if  any)  in  the  table  is  then  considered. 

The  printed  output  (Table  3)  f°r  each  mode  consists  of 
the  path  parameters,  the  reflections  that  take  place  for  the 
mode  in  "decoded”  form,  the  frequency,  take-off  angle, 
ground  range,  and  delay  tljae  for  the  rode#  ii*t  of  heights 
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vs  range  for  the  reflection  points  on  the  path  is  printed 
if  this  information  is  requested  of  the  program  (Table  4). 

This  process  continues  for  each  frequency  in  the 
specified  range  until  the  frequency  range  is  exhausted  or 
until  some  frequency  in  the  range  fails  to  propagate  any  rays 
between  the  transmitter  and  receiver.  At  this  point  new 
ionospheric  and/or  new  path  data  may  be  read  into  the  program 
and  the  process  may  be  repeated. 


B.  PROGRAM  DETAILS 

The  ray-tracing  program  is  written  in  the  FORTRAN  II 
computer  language,  specifically  for  an  IBM  '090  data- 
processing  system.  With  modifications  to  the  input  output 
statements,  the  program  could  probably  be  adapted,  with 
little  trouble,  to  other  systems  such  as  the  IBM  709,  CDC 
l604,  etc.  The  program  requires  two  magnetic-tape  units, 
one  designated  logical  unit  11  and  the  other  designated  as 
logical  6.  The  tape  on  logical  unit  11  serves  as  an  input 
tape  and  the  tape  on  logical  unit  o  is  trie  output  tape. 

Data  for  the  raj.  -tracing  program  are  prepared  by  a 
second  program,  which  shall  be  referred  to  as  the  data 
program".  The  lata  program;  requires,  as  Its  input,  atles 
of  f  Fg  and  M  3000  factor  coefficients  for  the  month  for 
which  ray  tracing  is  to  be  done,  in  adult Ion  to  parameters 
associated  with  the  paths  which  are  to  be  ray-traced.  The 
output  of  the  data  program  Is  a  magnetic  tape  tha4  Is  used 
as  the  Input  tape  for  the  ray-tracing  program. 

The  tables  of  f  F,  and  M  3000  factor  coefficients  (Dsk^ 
may  be  obtained  in  punched-card  form  from  the  Bureau  of 
Standards  CRPL  at  Boulder,  Colorado.  (See  CRPL  Ionospheric 
predictions,  Handbook  90.  [Ref.  1 8 ] )  These  tables  also  con¬ 
tain  the  sunspot  number  for  the  month. 
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The  data  program  will  also  accept  as  input,  actual 
measured  values  of  fQFg-  ond  F^-layer  real  height  a~  a 
function  of  their  position  on  the  path,  to  be  used  by  the 
ray-tracing  program. 

Appendix  F  contains  listings  of  the  FORTRAN  source  pro¬ 
grams  for  both  the  data  and  ray-tracing  programs  in  addition 
to  sample  output  from  the  ray-tracing  program  which  has 
already  appeared  as  Tables  2,3 ,  and  4,  and  sample  input  for 
the  data  program. 

C.  OPTIONS  AVAILABLE  ON  THE  DATA  PROGRAM 

In  preparing  data  for  paths  to  be  ray  traced  for  a  given 
month,  the  tables  (on  cards)  of  D=k  need  be  read  only  once 
by  the  program.  When  the  program  reads  the  Dsk  tables  from 
cards  it  places  them  on  magnetic  tape  (logical  unit  10),  in 
binary  form,  where  they  are  available  for  future  use. 

Normally,  ray  tracing  is  done  over  a  given  path,  for  a 
24-hour  period,  once  each  hour:  however,  provisions  have 
been  made  in  the  data  program  to  allow  tracing  for  an  arbi¬ 
trary  number  of  selected  times  (at  most  100)  on  any  given 
path.  The  number  of  times  to  be  traced  is  specified  to  the 
program,  followed  by  the  actual  times  to  be  used.  For  exam¬ 
ple,  in  the  normal  case  £4  times  would  be  specified  followed 

by  each  hour  from  0  through  23. 

When  empirical  data  are  to  be  used  for  ray  tracing  a 
given  path,  that  is,  ionosonde  records  of  fQF2-  and  F2' 
layer  real  heights,  even  though  the  CRPL  tables  of  Dsk 
are  not  used,  the  Dgk  tape  must  be  mounted  on  unit  10  none¬ 
theless.  These  empirical  data  are  presented  to  the  data 
program  in  the  form  of  f  Fg-  anc*  Fg-layer  height  as 
a  function  of  distance  along  the  great-circle  path,  measured 
in  kilometers  from  the  transmitter,  Ther»  must  be  an  odd 
number  of  measurement  points  specified. 


£7 
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Sporadic-E  data  may  also  be  included  when  the  empirical 
data  form  is  used.  No  provisions  have  been  made  uo  include 
sporadic-E  when  the  prediction  tables  are  used.  However, 
this  omission  may  be  remedied  with  only  minor  difficulty; 
and  procedure  will  be  discussed  after  a  description  of  the 
output  from  the  ray-tracing  program. 

The  data  program  writes  a  BCD  tape  on  unit  11  which  is 
used  as  an  input  tape  by  the  ray-tracing  program. 


D.  INFUT  TO  AND  OUTPUT  FROM  THE  RAY-TRACING  PROGRAM 


Input  to  the  ray-tracing  program  is  provided  by  a 
binary-coded  decimal  (BCD)  tape  written  by  the  data  program. 
It  is  to  be  mounted  on  unit  11.  Output  from  the  ray-tracing 
program  consists  of  an  ionospheric  profile  (i.atle  2),  con¬ 
structed  from  either  the  CRPL  predictions  or  empirical  data, 
at  100 -km  intervals  along  the  path;  path  identification 
information  such  as  the  name  or  number'  of  the  pa  >h ,  the  tirv  , 
month  and  day  for  which  the  tracing  is  being  done,  etc . 

The  actual  path-identification  information  used  is  p  to  th* 
user  and  will  t ■  explained  in.  the  section  or.  the  pi1'  paratior 


of  input  cards  for  the  data  program. 

The  length  of  the  great -circle  path  se tween  the  trans¬ 
mitter  and  receiver,  the  latitude  and  longitude  of  the 
transmitter  (Tx ) ,  the  bearing  from  the  transmitter  to  the 
receiver  are  all  printed  and  labeled  fox*  each  time  a  series 
of  rays  is  traced.  (In  the  normal  case,  once  each  hour  for 
the  24-hour  period.) 

Actual  information  concerning  the  rays  traced  appear^ 
in  Table  3  under  the  following  column  headings,  with  the 


associated  definitions: 

MODE;  The  mode  structure  of  the  ray  propagated  between 
the  transmitter  and  receiver.  The  symbol  -E 
Indicates  a  ray  reflection  from  the  bot tom  side 
of  the  Z  layer.  The  symbol  Indicates  a  ray 
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reflection  from  the  toD  side  of  the  E  layer.  The 
same  definitions  apply  to  F s  F,,  and  F2  layers. 
Obviously.  "-F2"  is  not  defined  and  will  not  occur. 

FREQ:  The  frequency  (in  megacycles)  of  the  ray  traced. 

BETA:  The  take-off  angle  (in  degrees)  of  the  ray  traced. 

DIST:  The  actual  ground  distance  (in  kilometers)  the  ray 

- *  travels  between  transmitter  and  receiver.  Because 

of  inaccuracies  in  the  ray-tracing  technique  this 
distance  will,  in  general,  not  be  equal  to  the  actual 
path  length. 

TIME:  The  delay  time  (in  milliseconds)  of  the  ray  traced, 
corrected  to  the  ac tual  path  length. 

DIFF:  The  difference  (in  kilometers)  between  the  ground 

distance  the  ray  travels  and  the  actual  path  length. 

DB:  The  attenuation,  (in  decibels)  the  ray  experiences 

'  due  to  D  layer  absorption  only . 

Information  concerning  reflection  heights  of  the  rays  is 
also  printed  (Table  4),  but  only  if  it  is  specifically  re¬ 
quested  of  the  program.  The  details  for  obtaining  this 
information  are  discussed  in  the  next  section. 

If  the  reflection-height  information  is  requested  it 
appears  in  the  following  form:  the  path  parameters  and 
identification  are  printed,  the  MODE  is  specified  (as  above), 
along  with  the  ray  frequency  and  take-off  angle.  The  heights 
appear  under  a  column  headed  HEIGHT  and  the  corresponding 
range  appears  under  a  column  headed  RANGE ;  both  are  in 
kilometers . 

In  all  cases  an  ionospheric  profile  along  the  path  is 
printed  prior  to  the  printing  of  any  other  information.  If 
consists  of  the  path-identification  information  and  columns 
headed  FOE,  FQF1,  F0F2,  HT  F0F2  and  RANGE.  The  fQ  values  are 
in  megacycles  and  the  fQF2  height  column  is  In  kHorrv  ter?  , 
as  are  the  ranges.  The  range  is  measured  from  the  trans¬ 
mitter  end  of  the  path,  and  the  values  fall  on  the  great 
circle  between  the  transmitter  and  receiver.  The  path 
parameters  printed  consist  of  the  PATH  LENGTH,  IX  LAT 
(transmitter  latitude),  TX  LONG  (transmitter longitude) , 
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and  RX  BEARING  (the  great-circle  bearing  from  the  transmitter 
to  the  receiver) . 

In  order  to  provide  sporadic -E  information  on  ray 
tracings  that  make  use  of  the  CRPL  prediction  tables,  it  is 
necessary,  first,  to  accomplish  the  required  tracings  with¬ 
out  Eq  data,  and  then,  using  the  fQP2-  and  Fg-layer  real- 

u 

height  information  provided  by  the  ray-tracing  program, 
resubmit  this  information  to  the  data  program  in  the 
empirical-data  format,  along  with  the  required  Es  data. 

This  technique  was  adopted  in  the  interest  of  programming 
simplicity  and,  since  the  inclusion  of  Es  is  usually  done  or* 
an  "after-the-fact"  basis,  it  would  seem  a  justifiable 
approach . 


E.  INPUT -CARD  FORMATS  FOR  THE  DATA  PROGRAM 


1 •  First  Card 

The  first  card  of  every  data  set  contains  the  pro¬ 
gram  variables  called  IDATA,  IDSKC,  IEND,  in  that  order. 
This  card  is  read  under  a  FORTRAN  format  of  (312).  The 
value  of  each  of  these  variables  may  be  "1"  or  "0"  (zero'). 

Tf  that A  =  1:  Data  are  to  be  supplied  to  the  pro¬ 


gram  in  the  empirical  format. 

If  IDATA  =  0:  Data  are  uo  be  supplied  to  the  pro¬ 
gram  In  the  form  to  mak-'  use  of  the  CRPL  D,jk  tables. 

If  IDSKC  =  1:  The  D^k  tables  for  the  month  in 
question  have  not  yet  been  put  on  magnetic  tape  and 

immediately  follow  this  first  card. 

If  IDSKC  =  0:  The  D3k  tables  for  the  month  in 

question  are  on  magnetic-tape  unit  10. 

If  IEND  =1:  An  END  OF  FILE  mark  is  to  be  written 
immediately  on  magnetic-tape  unit  11  and  program  execution 
is  to  be  terminated. 
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If  I END  •  Of  Additional  sets  of  path  data  follow. 


2.  Cards  for  Program  Using  -isk  Tables 

The  following  cards  constitute  the  information 
required  to  generate  data  for  the  ray-tracing  program  using 


the  Do^  tables: 


Card  #1:  Contains  the  program  variables. TXLAT,  TXLON, 
RXU'T’/'mON,  SUNDEC. 

TXLAT:  The  latitude  of  the  transmitting  point  in 
degrees  and  hundredths  of  degrees.  North  latitude 
is  +;  South  latitude  is  Format  F7.2. 

TXLON:  The  longitude  of  the  transmitting  point  in 

degrees  and  hundredths  of  degrees.  East  longitude  is 
+;  West  longitude  is  Format  F8.2. 

RXLAT :  The  latitude  of  the  receiving  point.  Same  as 
TXLAT . 

RXLON :  The  longitude  of  the  receiving  point.  Same  as 
TXLON. 

SUNDEC:  The  declination  of  the  sun  in  degrees, 

Format  F?.2. 

This  information  is  obtained  from  the  Nautical 
Almanac . 

Card  #2:  Contains  the  program  variables  FREQL,  FREQD, 

FREQH,  AU&LL,  ANGLD,  ANGLH. 

FREQL:  The  lowest  frequency  to  be  traced,  in  mega¬ 
cycles.  Format  FJ.3. 

FREQD:  The  frequency  increment  to  be  used  between  the 
lowest  frequency  and  highest  frequency,  in  megacycles. 

Format  F7.3.  . 

FREQH:  The  highest  frequency  to  be  traced,  m  mega¬ 
cycles.  Format  F7.3. 

ANGLL:  The  lowest  take-off  angle  to  be  traced,  in 
degrees.  Format  F7.3. 

ANGLD:  The  take-off  angle  increment  to  be  used  between 
the  lowest  angle  and  highest  angle,  in  degrees. 

Format  F7.3*  ,  ^  , 

ANGLH:  The  highest  take-off  angle  to  be  traced,  m 

degrees.  Format  F7.3 

Card  #3:  Contains  the  program  variables  NTIMES . 

- HTTHES:  The  number  of  specific  times  of  day  to  be 

used  by  the  ray-tracing  program  for  the  path  described. 
0  <  NTIMES  <  100.  Format  I  3. 

Card  Set  #4:  Contains  the  program  variable  TIME  (I).  One 
card  for  each  value  of  TIME  (),  in  GMT,  to  be  used.  The 
number  of  cards  must  correspond  to  NTIMES.  Format  Fb.2. 

Card  #5:  Contains  the  program  variable  NCHT .  NCHT  is 
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the  number  of  discrete  frequencies  for  which  detailed  ray- 
re  flecTTon-  height  information  is  desired.  If  no  such 
information  is  desired,  NCHT  =  0.  Format  1 4.  0  <  NCHT  <  50. 

Card  Set  #6:  Contains  the  program  variable  HFREQ,  (I), 

I  =  1 ,  2, . .  .“NCHS.  There  are  NCHT  cards  in  this  set,  each 
containing  a  discrete  frequency  for  which  ray- ref  lection- 
height- information  is  required.  If  NCHT  =  0  there  are  no 
cards  in  this  set.  Note  -  frequencies  specified  must  corre¬ 
spond  to  frequencies  specified  to  be  traced  by  the  progra:  . 
Format  F7.3. 

Card  Set  #7:  Contains  alpha-numeric  data  in  columns  1 
thru  60  i'orTdentificatlon  purposes.  One  card  must  appear 
for  each  time  used.  The  actual  information  used  is  at  the 
user's  discretion. 

Additional  sets  of  path  data  may  follow,  providing 
each  set  is  prefaced  by  a  card  as  described  in  Sec.  1. 

A  card  of  the  type  described  in  Sec.  1  with  IEND  =  1 

should  immediately  follow  the  last  set  of  path  data. 

3.  Cards  for  Program  Using  R^al -Height  Measurements 

The  following  cards  constitute  the  information 
required  to  generate  data  for  the  ray-tracing  program  using 
actual  measurements  of  fQFg-  anc*  F^-layer  real  height. 
Remember,  it  is  necessary  to  have  the  tape  mounted  on 

unit  10,  even  though  the  tables  are  not  used! 

A  card  as  described  in  Sec.  1  with  IIA.A  =  1. 

Card  #1:  Same  as  card  if  1,  Sec.  1. 

Card  fU :  Same  as  card  #2,  Sec.  1. 

Card  #3:  Contains  the  program  variable  NSETS .  N'SET.  : 

The  number  of  sets  of  empirical  data,  for  the  path 
described  by  cards  #1  and  #2,  to  be  read.  0  <  NSETC  . 

Format  I  3 . 

Card  #4 :  Same  as  Card  // (  ot  Sec.  1. 

Card  #5:  Contains  the  program  variables  NPTS,  SSN,  H  TJR. 
NPTS  *.  The  number  of  measurements  along  the  path  a^ 
described  by  cards  $1  and  #2. 

3  <  NPTS  <100  and  must  be  odd .  Format  13. 

SSN:  Sunspot  number.  Format  F5.1- 
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Forma  t 


HOUR:  The  time,  in  GMT ,  of  the  measurements. 

F6.2. 

Card  Set  #6:  Contains  the  program  variables  AF0F2(I), 

AHT(I),  DIST(I).  I  =  1,  2,  ...  .  NPTS. 

AF0F2(l):  f  Fp  in  megacycles  at  the  point  I.  Format 

F6.2.  °  ,  , , 

AHT(I):  Real  height  of  the  layer  maximum  at  the 

point  I,  in  kilometers.  Format  F7.2. 

DIST(I):  The  distance  from  the  transmitter  along  the 
great-circle  path  to  the  point  I,  in  kilometers.  Thf. 
first  measurement  must  be  at  the  transmitter 
(DIST  (l)  =  0)  and  the  last  measurement  must  be  at 
the"  receiver  (DIST  (NPTS)  =  path  length).  Format  F9-2 

Card  #7:  Contains  the  program  variable  IES. 

—  TfTES  =  0,  no  Es  data  are  to  be  considered. 

If  IES  =  1,  E  data  immediately  follow.  Format  13. 
s 


4.  Cards  for  Program  Using  Sporadic  E  Data 

Card  #1:  Contains  the  program  variable  N PATCH . 

— N PATCH •  The  number  of  sporadic  E  patches  on  the  path 
this  particular  time.  0  <  NPATCH  <  10.  Format  14. 

Card  Set  #2:  Contains  the  program  variables  PSTART(I), 

- PEND(DT  1  =  1*  2, ...NPATCH. 

PST ART ( I ) :  The  distance  from  the  transmitter,  along  the 
great-circle  path,  of  the  starting  point  of  E^  patch 
number  I. 

PEND  (I):  The  distance  from  the  transmitter,  along  he 
great- circle  path,  of  the  ending  point  of  E^  patch 
number  I.  Format  IX,  2F9.2. 

Card  #3 

fTpTT 

"I".  0  <  NPT  <  10.  Format  I  4. 


Contains  the  program  variable  NPT. 

The  number  of  f  E  values  to  be  read  in  for 
o  s 

o  < 


patch  number  "I".  0  <  NPT  _ 

Card  Set  #4:  Contains  the  program  variables  ESDIST  (I,J), 

TFOES  ( I, J)  1=1,  2. ...NPATCH,  J  =  1,  2,...NPT. 

FSDTST  (i.j)*  The  distance  from  the  transmitter  along 
the  great-circle  path  to  the  point  (I,J).  Note  that 
ESDIST  (1,1 )  must=FSTART ( I )  and  ESDIST ( I, NPT )  must 

M 

f _E„  at  the  point  (I,J).  Format  IX,  2F9.2 


=  PEND 
TFOES ( I, J): 


o  S 
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F.  DATA-SET  EXAMPLES 

Three  sets  of  examples  are  Included  at  the  end  of  this 
report  to  illustrate  graphically  the  preparation  of  data 
sets  in  the  form  described  above  In  Sec.  E. 

G .  PROGRAMS 

Card  decks  of  the  FORTRAN  source  programs  for  both  the 
data  and  ray-tracing  programs  are  available  from  the  Stanford 
Radioscience  Laboratory. 

H.  RUNNING  THE  PROGRAMS 

To  run  the  data  program,  prepare  the  appropriate  data- 

input  deck  in  the  appropriate  format,  and  submit  this  with 

the  7090  binary  deck  for  the  data  program,  along  with  the 

appropriate  control  cards  for  the  FORTRAN  MONITOR  in  use. 

(This  varies  with  the  7090  installation.)  Specify  the 

tapes  to  be  mounted  on  logical  units  1  and  11.  Naturally, 

if  D  tables  are  to  be  read  from  tape,  a  specific  tape  must 
sk 

be  mounted  on  unit  10.  At  the  end  of  the  run,  unit  11  will 
contain  the  input  data  to  be  used  by  the  ray-tracing  pro¬ 
gram. 

To  run  the  ray-tracing  program,  mount  the  appropriate 
data  tape  on  logical  unit  11.  Submit  the  7090  binary  deck 
for  the  ray-tracing  program  along  with  the  appropriate 
FORTRAN  MONITOR  control  cards.  Output  from  this  program 
appears  on  the  "normal"  FORTRAN  output  tape  unit  f  . 

Neither  program  makes  use  of  any  sense  switches  or  other 

console  features. 
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VII .  CONCLUSIONS 


In  describing  the  reasons  why  the  Kift-Fooks  technique 
was  chosen,  how  it  would  be  used  in  the  analysis  of  prop¬ 
agation  data,  and  giving  details  of  the  program  for  use  on 
a  high-speed  digital  computer,  no  comparisons  were  made  with 
actual  records  taken.  It  remained  the  intention  of  the 
authors  to  outline  the  work  done  here  at  Stanford  and  their 
reasons  for  doing  it. 

Comparison  of  ray  tracings  with  experimental  data  has 
been  done  on  several  paths  and  the  results  of  these  com¬ 
parisons  are  scheduled  for  another  report. 

Hopeiully,  the  reader  of  this  report  will  find  suffi¬ 
cient  information  to  enable  him  to  reproduce  this  version 
of  the  ray-tracing  program  for  use  on  available  computers 
should  he  desire  to  do  so.  Duplicate  decks  of  the  program 
can  be  obtained  from  the  authors  by  written  request. 
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APPENDIX  A.  TERMINOLOGY 


This  Appendix  has  been  taken  directly  from  the  Report 
of  the  Lindau  Meeting  on  Oblique  Sounding  of  the  Ionosphere,” 
May  6-10,  1963.  The  recommendations  listed  below  are  to  be 
submitted  to  URSI  and  are  also  scheduled  to  be  published  it. 
the  IQfY  notes. 

It  was  recognized  that,  for  purposes  of  data  interchange, 
a  need  exists  for  the  standardization  of  certain  terms.  As 
a  first  step  in  this  direction,  the  following  recommendations 


are  made . 


Canital  letters  should  be  used  in  oblique-incidence 
work  in  contrast  to  the  small  letters  agreed  upon 
in  vertical-incidence  work. 

In  view  of  the  ambiguity  in  the  meaning  of  «  usable  », 
the  term  maximum  usable  frequency  (MUF)  shoull  be 
eliminated  in  the  description  of  oblique-incidence 
ionograms . 

The  use  of  the  word  "virtual  path"  should  refer  to  the 
time  of  flight  (group  delay)  in  oblique  propagation 

work. 

In  ray  tracing  the  following  symbols  are  suggested 
( Fig  A1 ) 

V  a.’  0  *  the  angle  of  incidence  at  the  bottom  of  the 

0  ionosphere.  .  .  , 

b  ( 1)  the  angle  of  incidence,  at  the  real  height 

r  of  reflection,  of  the  extension  of  the  linear 
ray  path  below  the  ionosphere. 

c .  L  the  angle  between  the  ray  path  and  the 
vex-tical  at  any  point  along  the  path. 


A  Moo  5.rurlf  nf  PlC 
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The  following  terminology  is  suggested  for  the 

description  of  path  structure  (Fig.  A2). 

t  For  propagation  paths  involving  reflections  by 

different  layers,  the  reflections  (or  hops)  should  be 
specified  in  order  of  their  position  with  respect  to 

the  transmitter.  Thus  bE  -  3F2  ^^^Lsmitter 
reflections  from  the  E  layer  near  the  transmitter 

followed  by  three  reflections  from  the  F^  1  y  ( 

A2-a ) . 


F2F2 


FIG.  A2  RECOMMENDED  MODE  IDENTIFICATION. 


-  39  - 


^el-63~103 


6.  The  use  of  a  dash  is  convenient  for  the  representation 
of  a  ground  reflection.  The  absence  of  a  dash^will 
then  show  up  M-type  ray  paths  and  1  supermodes, 
example  F  -  Es  -  F2  (Fig.  A2-b)  represents  an  F-layer 

hop  followed  by  ground  reflection  to  the  lower  side 
of  the  E„  layer,  reflection  back  to  ground,  then 

reflection  to  the  lower  side  of  the  F2  layer  and  finally 
back  to  ground.  On  the  other  hand  F  Eg  F2  (Fig.  A2- 

represents  an  M-type  path  in  which  the  ray  is  reflected 
from  the  F  layer  to  the  upper  side  of  the  E-.  layer, 

back  up  to  the  lower  side  of  the  F2  layer  and  down  to 
the  ground.  The  symbol  F2F2  (Fig.  A2-d)  means  an  F2 
reflection  followed  by  another  F2  reflection  without 
an  intermediate  ground  reflection  (supermode). 


The  following  terms  are  suggested  for  the  description 

of  oblique  ionograms  (Fig.  '3). 

7.  MOF  (Maximum  Observed  frequency)  means  the  highest 
frequency  on  which  the  sounder-transmitter  signalo  arc 
observe!  on  the  ionogram,  regardless  of  the  propagation 
path  involved . 

8.  LOF  (Lowest  Observed  Frequency)  means . the  lowest 
frequency  on  which  the  sounder-transmitter  signals  are 
observed  on  the  ionogram,  regardless  of  the  propagat j.or. 
path  involved . 


These  terms  (MOF  and  LOF)  may  be  used  also  to  describe 
identifiable  modes.  For  example  2F2  LCF  means . the . 
lowest  frequency  (observed  on  the  ionogram)  which  is 
propagated  by  two  reflections  at  the  F2  layer  and  an 
intermediate  ground  reflection.  .The  2F2  MOF  mears  the 
highest  observed  frequency  associated  with  two-hop 
F2  propagation,  regardless  of  whether  the  signal  n 
propagated  by  refraction,  by  scatter,  or  by  a  com¬ 
bination  of  both  mechanisms. 

The  lowest  observed  frequency  of  the  high-angle  ray 
may  be  distinguished  from  that  of  the  low-angle  ray 
by  the  letters  H  and  L  respectively.  Thus  sF^  HI  F 
is  the  lowest  frequency  (observed  on  the  ionogram)  o 
the  signal  that  is  propagated  via  the  high-angle, 
two -hop,  F2  path  and  2F2  LLOF  is  the  lowest  frequency 
(observed  on  the  ionogram)  of  the  signal  that  is 
propagated  by  the  low-angle,  two-hop.  Fa  patn. 
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11. 


The  one -hop  modes  do  not  need  the  number  l(one)  in 
front.  For  example,  F2  LLOF  means  the  low-angle  ray 
LOF  for  the  one -hop,  F2  ray  path, 

12.  When  it  is  required  to  distinquish  between  the  ordinary 
and  extraordinary  ray  paths  an  "o'*  or  "x"  nay  follow 

in  parentheses.  The  F2  MOF(x)  is  the  maximum  observed 
frequency  of  the  extraordinary  wave  that  is  reflected 
once  at  the  F2  layer. 

13.  Often  the  MOF  for  an  identifiable  path  is  greater  than 
the  frequency  on  which  the  regularly  refracted  com¬ 
ponents  of  the  high-and  low-angle  rays  Join.  It  is 
suggested  that  the  latter  frequency  be  called  the 
"junction  frequency"  and  that  it  be  denoted  by  JF. 


FIG.  A3  RKCOMMKNDKI)  I ONOGRAM- SCAI. 1 NG  PAHAMKTKR.S . 
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APPENDIX  B.  CALCULATION  OS  THS  SUN»S  ZENITH  ANGLE,  x 


cos  X  =  sin  \  sln  \  *  003  \  OOS  ^2  C0S  (02  "  0l) 


cos  a 


sin  7^  =  sin  ? \q  cos  (^)  +  cos  \  3iri  (^)  cos  a 

cot  ^0^  -  0Q  j  =  ^sln  cos  a  -  cos  cot  sin  a 

e2  =  2tt  (neglecting  equation  of  time)  (g.i) 


where 

X  =  latitude  of  transmitter 
o 

0q  =  longitude  of  transmitter 

T  =  time  in  hours  (U.T.) 

A,  =  declination  of  sun 

a  =  bearing  S.  of  N.  of  receiver  from  transmitter 

Or,  =  longitude  of  sun 

"K  =  latitude  of  point  on  path 

0^  =  longitude  of  point  on  path 

d  =  distance  from  transmitter  to  point  on  path 
R  =  radius  of  earth 

The  data  program  computes  a  and  the  path  length,  using 
the  latitude  and  longitude  of  both  the  transmitting  and 
receiving  points  and  supplies  the  ray-tracing  program  with 
these  parameters,  In  addition  to  the  latitude  and  longitude 
of  the  transmitter.  The  data  program  also  supplies  a  set 
of  distances  d^^  at  roughly  every  500  km  along  the  path  at 
which  fQF2  and  F2  4000  MUF  are  supplied  by  the  data  program. 
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APPENDIX  C.  A  METHOD  FOR  COMPUTING  F2  LAY£R  HEIGHT  h  FROM 
VALUES  OF  fQF2  AND  F2  4000  MUF  m 


A  r.omogram  of  height  hm  versus  the  ratio  of  F2  4000  MUF 

and  f  F2  for  a  parabolic  layer  with  Y  =  0.4  h  is  presented 
o  ra  o  * 

in  the  Fooks  report  [Ref.  12],  A  polynomial  expression, 

valid  for 


2.15  < 


F2  4000  MUF 

- TT2 - 

o 


<  4.09, 


(C.l) 


which  approximates  the  nomogram  with  maximum  error  in 
h^  of  ±  6  km,  is  used  in  the  program  to  compute  h  . 

Let  x  =  ~--xr°T?  *  1-1  (^3000), 

o 

hm  =  ( ~TT7ui§)+  19,44  (4,09  -  x)  (x  "  2*15)  ( 

'x  '  +  46.0  (3.0  -  x)  (4.09  -  x)  (x  -  2.15) 
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APPENDIX  D.  CALCULATION  OF  REFLECTION  HEIGHTS  OF  THE  RAY  IN 

A  LAYER 


The  height  of  reflection  is  calculated  as  follows: 


h  =  h  +  Y 
r  o  m 


[>-f  Yf_gosij2J  (D  .1) 


in  the  case  of  a  ray  reflecting  from  the  bottom  of  a  layer, 
and 

f  /  \ 

(D.2) 


h  =  h  +  2Y  -  Y 
r  o  m  in 


[•■vnr 


cos  i 


in  the  case  of  a  ray  reflecting  from  the  top  of  a  layer. 

The  definition  of.  the  parameters  is  the  same  as  in 
Eqs.  (4)  and  (5). 


-  w  - 
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APPENDIX  E.  CALCULATION  OP  RAY  ATTENUATION  DUE  TO  D-LAYER 

ABSORPTION 

The  following  expression,  taken  from  RPU  No,  9 
[Ref.  21]  is  an  estimate  of  the  absorption  in  the  D  layer 

_  615.5  (1.0  +  0.0037-SSN)  .  cos1*3  (0.88lx) -N-sec  cpD 

DB " - (f  C^F98  (e.i) 

whe  re 

SSN  =  sunspot  number 
X  -  sun's  zenith  angle 

N  =  number  of  ray  passages  through  the  D  layer 
f  =  ray  frequency 

<PD  =  vertical  angle  which  ray  makes  with  the  D  layer 
fh  »  gyro -frequency 

DB  =  number  of  declbel3  of  ray  attenuation 

In  the  ray-tracing  program  the  D  layer  is  assumed  to 
be  at  a  height  of  70  km.  This  height  plus  ray  take-off 
angle  allows  the  calculation  of  q>D.  Since  the  program 

assumes  a  constant  ray  take-off  angle,  this  quantity  cpp 

is  computed  only  once  for  each  mode.  An  average  value  of 
cos  x  is  used  for  each  mode  and  an  average  value  of  f ^ 

along  the  path  is  used  in  the  calculation. 


APPENDIX  F  LISTING  OF  DATA  AND  RAY-TRACING  PROGRAMS,  SAMPLE 

OUTPUT  AND  INPUT  FORMATS 

The  following  figures  consist  of  sample  input  data  for 
the  data  program,  a  listing  of  the  data  program,  a  listing 
of  the  ray-trace  program,  and  sample  output  from  the  ray- 
trace  program. 
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6.0G  1.18  80)9.09  ? 7. 10  -6.26  669.57 

.E  6.00  13.19  8061.61  27.98  -'*.96  896.86 

5.00  1.22  8037.90  27.09  -7.65  617.28 

6. Or  1.28  8036.66  27.09  -8.69  290.56 

7.00  1.35  8035.68  27.09  -9.87  21  3.86 

7. GO  12.50  8026.50  27.86  -18.85  275.91 

8.00  1.66  60)6.67  27. C9  -10.78  163.91 

9. CO  1.56  8036.36  27.  10  -1  1  .01  129.56 

10.00  1.68  80)5.63  27.12  -9.92  1G6.J6 

1G .00  23.19  8015.76  30.17  -?9.59  66.07 

11.00  1.85  8039.26  27.16  -6.09  86.68 

il.00  13.20  7989.22  28.10  -56.13  60.36 
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25. GO  6.71  8016.11  27.91  -29,26  10.27 
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